function params_return= time_sampling(params)

dynatime = params.dynatime;
params_return = params;

if dynatime < 0
    error('dynatime must be larger than 0 !');
end

if dynatime ~= 0 % 时序采样
    D = params.D;
    T = params.dynatime;
    samplehz = params.samplehz;
    resol = params.resol;
    
    % sampling strategy
    sample_N = samplehz*T; % 总的采样点数
    delta_t = T/sample_N;
    t = (0:sample_N-1)*delta_t;
    
    % 只需要求半边功率谱，便于构造实数序列的双边功率谱
    % 所以只取一半时间频域
    t_miu = (0:(sample_N/2))/T; % time frequency
    f_y = (-resol/2:resol/2-1)/D; % spatial frequency
    
    % params_retun = params;
    params_return.t = t;
    params_return.t_miu = t_miu;
    params_return.f_y = f_y;
    params_return.sample_N = sample_N;
    params_return.delta_t = delta_t;
end

% end of the function
end